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We generally investigate necessary conditions for the generation of the long-distance entanglement 
of two spins indirectly interacting through a media system. First, we infinitesimally weaken the 
coupling between each focused spin and the media system in order to generate the long-distance 
entanglement. We give two necessary conditions for the media system to generate the long-distance 
entanglement in the weak coupling limit. We prove that the indirect interaction cannot generate 
the entanglement if it is 'classical.' We also give a necessary condition for the effective fields to 
satisfy for the generation of the long-distance entanglement. Second, we generate the long-distance 
entanglement by the use of only external fields. We show that external fields on adjacent spins are 
necessary in addition to external fields on the focused spins. A specific combination of the external 
fields creates the situation where the effective coupling between each focused spin and the media 
system is infinitesimal. Finally, we consider the cases where the coupling strength between each 
focused spin and the media system is finite. In particular, we show two examples where the external 
fields on the media system highly enhance the long-distance entanglement. 



I. INTRODUCTION 

The generation of the quantum entanglement [l[ is of great importance because the entanglement plays an essential 
role in quantum information processing 0, Many people have sought efficient generation of the entanglement 
between remote systems Q. When two systems, which we refer to as the focused systems, indirectly interact with 
each other via a media system (Fig.Q}, the entanglement of each focused system with the media system would decrease 
the purity of the focused systems. As a result, the entanglement between the focused systems rapidly decays and 
vanishes as the focused systems get separated distantly. Indeed, we cannot usually obtain the entanglement between 
systems far apart even in the low-temperature limit [f| @. For the generation of the entanglement between the 
remotely separated focused systems, we therefore need to suppress the entanglement between each focused system 
and the media system. It then appears to be a dilemma; we need to weaken the interaction between each focused 
system and the media system, and yet we need to keep the indirect interaction between the focused systems. 

Recent studies Q, however, reported successful generation of the entanglement between a far-separated spin pair 
in specific models; such an entanglement is referred to as the long-distance entanglement. In their models, the 
entanglement with the media system is suppressed enough although the two spins indirectly interact with each other 
via the media system. After the first paper on the long-distance entanglement, various systems have turned out to 
be usable for generating the long-distance entanglement In order to utilize the long-distance entanglement for 

practical applications, however, there are still many problems to solve. In the present paper, we find answers to the 
following questions on the long-distance entanglement: 

1. What are the conditions for a media system to generate the long-distance entanglement? 

2. What are the conditions to generate the long-distance entanglement by the use of only external fields on the 
focused systems? 

3. Is it possible for external fields on the media system to enhance the long-distance entanglement? 
These questions are closely related to practical applications. 
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FIG. 1: Indirectly interacting systems. The two focused systems interact with each other through the media system. The 
entanglement between each focused system and the media system decreases the purity of the focused systems. 
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FIG. 2: (a) A schematic illustration of the generation of the long-distance entanglement for Section Hill and IVl Two spins are 
connected to the media system with the coupling Hamiltonians Hi t media and -Hjv.media- We refer to the spins 1 and TV as the 
focused spins and the entanglement between the focused spins as the long-distance entanglement. In order to generate the 
entanglement between the focused spins, we decrease the amplitude of the coupling Hamiltonian. We consider infinitesimally 
weak couplings in Section iHll and finite couplings in Section IVl 

(b) A schematic picture of the generation of the long-distance entanglement by local fields for Secition IIVI We connect four 
spins to the media system and consider the entanglement generation between the focused spins 1 and JV. We apply the local 
fields to the four spins. By increasing the local fields h-z and ftjv_i, we can decrease the effective coupling between each of the 
focused spins and the media system. 



One of the most popular methods of the generation of the long-distance entanglement is to weaken infinitcsimaly 
the coupling between each focused system and the media system 04HI- It has been shown for several models that 
we obtain the maximum entanglement between the far-separated spins in the limit of the weak coupling. In some 
models, analytical calculations of the long-distance entanglement are possible and important quantities such as the 
first excitation energy are obtained in Refs. (8l-[l0j. In Ref. @, a mathematical analysis is also given as to why we can 
generate the strong entanglement in this framework. However, there is no analytical argument on a general condition 
to generate the long-distance entanglement, which is one of the most important problems in discussing the usefulness 
of the long-distance entanglement. 

It is also an important question whether we can generate the long-distance entanglement with external fields or not. 
It is a popular problem how much entanglement can be generated by modulating system parameters [5|, 1 1 81 - 4261 ] . The 
external field is one of the easiest parameters which we can control freely. Reference [2(| discussed the generation of 
large entanglement between distantly separated two spins in an XX spin chain; they showed that the local control of 
the external fields is an effective method of generating a large entanglement between the separated two sites in this 
specific spin chain. In general systems, however, we cannot generate the long-distance entanglement if we modulate 
the external fields only on the two focused spins. 

Answering the above questions, we discuss in the present paper the generation of the long-distance entanglement 
between the two focused spins which indirectly interact with each other through a media system. Our main results 
are summarized as follows: 

1. We first consider the system in Fig. [3] (a) and decrease the coupling strength between the media system and 
each of the focused spins 1 and TV. In the weak coupling limit, we give two necessary conditions for the media 
system to generate the long-distance entanglement in the forms of two sufficient conditions for the non-existence 
of the long-distance entanglement. First, the indirect interaction must not be 'classical.' Second, the effective 
fields on the focused spins must not be much larger than the indirect interaction between these two spins. We 
will give mathematical expressions of these conditions below. 

2. We next consider the generation of the long-distance entanglement with external fields. As has been expected, 
it is not enough to control the external fields only on the two focused spins. In addition to the external fields on 
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the focused spins, we also have to control the external fields on the system adjacent to the focused spins, as is 
depicted in Fig. [2] (b). We give a necessary condition for the external fields to give rise to infinitesimal effective 
couplings between the media system and each of the focused spins 1 and N, thus generating the long-distance 
entanglement between the focused spins owing to the mechanism given in the previous item. 

3. We finally consider the quantum system in Fig. [2] (a) with finite couplings. We can highly enhance the long- 
distance entanglement by modulating the external fields on the media system. This comes from the fact that 
the indirect interaction depends on the external fields. 

We show these results analytically and numerically. 

This paper consists of the following sections. In sectionHU we review the general framework of the generation of the 
long-distance entanglement. Sections ITlTl IIVI and \V\ describe the items 1, 2 and 3 above, respectively. In Section IVl) 
a discussion concludes the paper. 



II. GENERAL FRAMEWORK OF THE GENERATION OF THE LONG-DISTANCE ENTANGLEMENT 



In the present section, we review the general framework of the generation of the long-distance entanglement [7H11| . 
We consider a quantum system in which two 'focused' spins are connected to an arbitrary media system (Fig. [5] (a)); 
that is, 

-fftot = Hi n t + Hlf, (1) 

with 

H'mt = -//media ~\~ -//couple, 

-//couple = ^2 (fl (8 -ff^modia + a N ® ^/V.media) ) 

H LF = J2 (hi4 + h%4f) (2) 

where {a l }i= x , y ,z are the Pauli matrices of s = 1/2 spins, while {H{ media , H l N lnedia }i= x , y , z and //media are arbitrary 
(but not the identity) 2 JV ~ 2 -dimcnsional Hamiltonians which do not include the spins 1 and N. We assume that the 
ground state of H me dia is not degenerate. We hereafter refer to the spins 1 and TV as the focused spins and to the 
external fields {h\, h l N } i=x y z as the local fields. We also refer to the entanglement between the focused spins 1 and 
N as the long-distance entanglement. In fact, the long-distance entanglement is usually defined as the entanglement 
between infinitely separated two spins Q; if the distance between the two focused spins is large but finite, the 
entanglement between such a spin pair is often referred to as the quasi- long-distance entanglement [8| . In the present 
paper, however, we define all the entanglement which is generated by the indirect interaction as the long-distance 
entanglement for simplicity. 

Let us consider the problem of enhancing the ground-state entanglement between the spins 1 and N by modulating 
the coupling Hamiltonian H coup \ c . The density matrix of the total system in the ground state is given by 

e -/3fftot 

Ptot = lim , (3) 

P^oo Ztot(p) 

where /? is the inverse temperature (/cbT) -1 with fee the Boltzmann constant and Z to t(/3) = tr(e~P Htot ) is the partition 
function. The density matrix of the focused spins 1 and TV is 

piiv = triivptot, (4) 

where truv denotes the trace operation on the system except the focused spins 1 and N. In order to quantify the 
entanglement, we here adopt the concurrence |27| as an entanglement measure. The concurrence C(pin) is defined 
as follows: 

C(puv) max(Ai - A 2 - A 3 - A 4 , 0), (5) 
where {X n }f l=1 are the eigenvalues of the matrix 

\J PijvK ® cr y N )p* 1N (af <g) cr v N ) (6) 
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FIG. 3: Comparison of the entangled ground state and the non-entangled ground state. In the limit of weak coupling, the 
focused spins 1 and N cannot entangle strongly with the media system because of the monogamy. Therefore, we can consider 
the two cases (a) and (b). In the case (a), the energy of the entangled state is less than the non-entangled state, while in the 
case (b), the energy of the non-entangled state is less than the entangled state. 



in the non-ascending order Xi > X 2 > A3 > A4. 

A popular method of generating the long-distance entanglement in the ground state is to reduce the coupling 
Hamiltonian Hcoupie to zero [7HTlj|: namely 



where j| ■ • • | 
limit © as 



ll^couplol! ->■ 0, (7) 

denotes the matrix norm. It has been shown that the long-distance entanglement is generated in the 



lim C{puf) ? 

Uncouple |j— >0 



(8) 



in some quantum systems [7H11|. for example, an XX spin chain. 

Let us explain qualitatively why we can generate the long-distance entanglement in the limit ([?]). A key feature 
for the generation of the long-distance entanglement is the monogamy [28j . which is a property that multiple pairs 
cannot share a strong entanglement simultaneously. Because of the monogamy, if the components of the media system 
are entangled strongly with each other, a large entanglement cannot be generated between each focused spin and the 
media system. In the limit ([7]) of weak coupling, the components of the media system may be entangled strongly with 
each other, but are scarcely entangled with each focused spin. Then, the focused spins can be entangled strongly only 
with each other. 

It then appears to be claimed that the focused spins are entangled without the interaction, which is obviously not 
the case. The fact is that the interaction is weak but not zero, and hence the quantum correlation between the two 
focused spins can decrease the energy of the total system. Then the ground state may support the entanglement 
between the focused spins (Fig. [3|. We can describe the indirect interaction in terms of the effective Hamiltonian 
in the limit ([7]) . We will give the effective Hamiltonian in Section MI Bl in order to obtain general properties of the 
long-distance entanglement. 

The next question is how to realize the weak-coupling limit (J7J in realistic situations. We answer it by modulating 
local fields in the system in Fig.[2](b), extending Fig. [2] (a) slightly. By increasing the local fields on the spins 2 and 
N — 1, we can effectively decrease the coupling between the media system and each of the focused spins 1 and N. In 
fact, as shown in Section llVl we cannot achieve the long-distance entanglement by simply letting \h 2 \ and |/ijv-i| to 
infinity. The merit of using the external fields is that we can more easily and precisely control them than the coupling 
interactions. We consider the following Hamiltonian: 



= E ( 



media 3 



(9) 



i—x,y,z 



where the Hamiltonian -ff m cdia is now an arbitrary 2 -dimensional Hamiltonian and we assume that the ground 
state of -ffmcdia is not degenerate. We modulate only the local fields {h\, h\, h l N _ 11 h l N } i=x ^ y ^ in order to generate the 
entanglement between the spins 1 and N . In Section ITVl we clarify the exact correspondence between the decrease of 
the interaction and the increase of the local fields. 
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III. GENERAL CONDITIONS FOR THE GENERATION OF THE LONG-DISTANCE 

ENTANGLEMENT 

In the previous section, we presented general statements on the generation of the long-distance entanglement. 
However, they do not tell us whether the media system contributes to the generation of the long-distance entanglement 
or not. We here present two cases in which wc cannot generate the long-distance entanglement. In the first case, 
the indirect interaction cannot generate the entanglement at all; we refer to such an interaction as the 'classical' 
interaction. In the second case, the effective fields on the focused spins seriously destroy the entanglement. These 
constitute sufficient conditions for the non-existence of the long-distance entanglement, thereby giving two necessary 
conditions for its existence as their contrapositions. Throughout the present section, we consider the system fl} in 
Fig. 13 (a). 

A. Classical interaction 

The first necessary condition for the generation of the long-distance entanglement is that the indirect interaction 
between the spins 1 and N must not be a 'classical' interaction. We define that the indirect interaction between the 
focused spins 1 and N is 'classical' if there exists the following separation of the Hamiltonian H lnt '- 

H in t = H A (cr 1 ) + H B (cr N ) 

with [H A {a 1 ),H B {a N )} = 0. (10) 

Note that Ha does not contain nor Hb contains o\. We can prove the following theorem. 

Theorem 1: If the interaction Hamiltonian H ln t is 'classical,' we cannot generate the entanglement between the 
spins <j\ and o"jv for any values of the local fields {h\}i= Xt y iZ and {h l N }i =Xt y tZ ; namely 

C( PlN ) = 0, (11) 

where the density matrix piN is defined in Eqs. ([3]) and Q and the concurrence C(puv) is defined in Eq. ([5]). 
Comments: For example, we can separate the following Hamiltonian in the form (|10D : 

TT TZ Z Z I JX X X I JX X X { 1 c~t A 

H int ^ J 1 o- 1 a 2 + J 2 a 2 a 3 + J 3 a 3 a 4 , (12) 

where the spin pair (1,2) interact with each other through the Ising interaction along the z-axis, while the spin 
pairs (2, 3) and (3, 4) interact with each other through the Ising interaction along the rc-axis. We can separate this 
Hamiltonian into Ha{<?i) and H 5(04) as 

H B (a 4 ) = JZofot. (13) 

These Hamiltonians satisfy the condition [Ha(o~i), Hsicr^)] = 0, and hence we cannot generate the entanglement 
between the focused spins <7i and 174 in this system for any values of the local fields {h\}i =x ^ y ^ z and {h\}i =x . y z . Note 
that the spins 1 and 4 are classically correlated with each other. If wc replace Jf erf erf in Eq. (|T2|) by Jf erf erf as 

Hint = Jl°~lV2 + ^ G 2 G l + ^3 CT 3 cr 4i ( 14 ) 

we cannot separate the Hamiltonian H lnt into the forms of Ha{o~i) and Hb{o~4) which satisfy [Ha{o~i) 1 Hb^^)] = 
anymore, and hence the spins 1 and 4 can entangle with each other. It is worth noting that the coupling in the 
z direction to the media spin 3 makes the 'classical' interaction (|13p into the non-classical one. 

Eigenstates can have the entanglement even if the condition (| 10[) is satisfied. For example, the Hamiltonian for the 
focused spins 1 and 3, 

-fftot = Hint + h*Oi + h 3 a 3 , 

Hint = J*0\(*2 + ^2 a 2 a ?, (15) 

with h\ = /if = Jf = Jf satisfies the condition (fTU|) . but it has an eigenstate (IT1T2T3) + ll^ls))/"^, which is highly 
entangled. Mixing of all the eigenstates with the Boltzmann weight always destroys the entanglement between the 
focused spins in the case of Eq. (|15|) . 

Finally, under this condition, for appropriate values of the local fields wc can generate the quantum discord, which 
is one of the non-classical correlations. Wc discuss the quantum discord in Appendix [a1 
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Proof: We prove the following equality under the condition fj 10[) : 

tn N (e-^)=Y^p^p%, (16) 

where truv denotes the trace operation on the system except the focused spins 1 and N , and the density matrices p™ 
and p N are physical states, namely, positive matrices. Then, the spins 1 and N are not entangled with each other by 
definition. By proving Eq. (|16[) , we can prove in the limit /3 — > oo that the density matrix (j4|) is decomposed into the 
mixture of the product states. 

First, under the condition (jlOp . we can decompose the density matrix as follows: 

e -0H tot _ e -/3HA(<ri) e -/3H B (<Tiv)^ Q 7 \ 

We can express e -0 H A(a) an( j e -PH B (a N ) ag 

-PH A (cn) _ V - ^ a „ In „ r 
e — / „ °1 <& Pmcdia <» -»iV, 

^t— 0,x,y,z 

e -PH B (*„ )= /l®/^®^, (18) 

i^— 0,x,y,z 

where 7i and In are the identity matrices in the spaces of the spins 1 and N, respectively, and we define a® = I\ 
and a N = In- We also define that the matrices p*f odia and p^dia are Hermitian operators in the media-spin space. 
Because Ha{ti) and Hb{cn) commute with each other, e~@ HA ( Ul ) and e~^ HB ^ N > also commute with each other. 
Therefore, we obtain the following equation: 

tr 1JV [a? <g> a ^ e -mA{^) e -m B {a N )^ = tT iN ^ g ^ e -/3H B ( CTN ) e -^. 4 ( CT1 )] , (19) 
where tr 1Ar denotes the trace operation only on the spins 1 and N. From this equation we obtain 

rmediarmedia — "media"media 

for p,v = 0,x,y,z. Therefore, the matrices pj^^ an< ^ Pmedia have simultaneous eigenstates. Then, we can express 

Pnfcdia an d Pmedia aS 



and 



/^mcdi; 



Nu 
"media 



X A^|n,^i,i/j V )(n-,Mi,i / A r l ( 21 ) 



= X r"|Ti,^i,i/jv)(n,/ii,i'jv|, (22) 

n=l 



where {|n, /ii, fjv)} are 2 W 2 pieces of the simultaneous eigenstates of pj^ edia and Pmedia- As a result, we obtai 



e -j3H A { ai ) e -l3H B {a N ) 



= £ A m t >i® KMi.^X^i."*!®^. (23) 



where the indices // and f' in the first line can be arbitrarily chosen (i/, /i' = 0, x, y, z), and hence we choose p! and 
z/ in accordance with p and za By tracing out the media spins, we have 

= E(E*X)«(E ( 24 ) 

n n v 
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Do not commute 

-(1 Commut^^Commute i)- , 

/^media Pmedia r media 
\n ,fii,p N ) — : — : — : — \no,n[,u N ) 

\ni,/j,i,v N ) — '- — — '■ — '■ — \ni,fj>' x ,u N ) 

\ |ni + 1,,%, i/jv) 



|m +2,/ii,i/ iV > ^ |m + 2,^;,^) 

|ni + 3, Hi,u N ) |m + 3,/iiji/jv) 

FIG. 4: A schematic picture of the eigenstates of pj^ edi& , Pmcdia an d P^dia- Suppose that the noth and the nith eigenvalues of 
Pmcdia are n °t degenerate (r"° 7^ r" 1 ); then the eigenstates |no, /xi , i/jv) and |m, /Ui, i*jv) are orthogonal to each other. Suppose 
that all the eigenstates of p^odia are degenerate at the (m+l)th and the (ni+2)th levels, namely, r™ 1 + 1 = r" 1+2 for v — 0,x, y, z. 
Then, any superposition of the states \m + 1, pi, Vn) and |m + 2, /Lti, i/jv) can be the (m + l)th eigenstates of Pmedia, an d hence 
it is possible that p^ edia and p^ dia do not have the simultaneous eigenstates in the space {|ni + 1, /ii,i/jv), |tii + 2, /zi, i'jv)}- 
We then form a block which is composed of the states |m + 1, pi, vn) and \m + 2, /zi, i^jv). 



At this moment, we cannot say that A"c>i and ^ r v™o v N are necessarily physical states, namely, positive matrices. 
In the following, we prove that Eq. (|24p can be reduced to the mixture of the product states as in the form ([TB")) . 

For the purpose, we should pay attention to the degeneracies of the matrices p^, dia and Pmcdia- ^ n f ac *> if there 
are no degeneracies in the eigenspaces of all these matrices for p, v = 0,x,y,z, we can easily prove that each of 
an d J2v T v a N ( n — 1,2,..., 2 JV_2 ) in Eq. (f2~4"|) is a positive matrix. Since the matrices p]^ cdia and p^ dia 

commute with each other as well as p^ cdia and p^ dia do, the matrices p^, dia and p^ cdia should also have simultaneous 
eigenstates if there are no degeneracies. If there are absolutely no degeneracies in all eigenspaces of p]^ cdia and pmcdia 
(p, v = 0, x, y, z), we have an orthonormal set of 2 N ~ 2 pieces of states |n), each of which is the simultaneous eigenstatc 
\n, p\, vn) for all of p, v = 0, x, y, z. Then, we have from Eq. (|18|) 

J2 A«<=tr liV (e-^^)|n)(n|), (25) 

fi— 0,x : y,z 

J2 Ty N =tr 1N (e-? H "^\n)(n\), (26) 

v— Q,x,y,z 

for n = 1, 2, . . . , 2 N ~ 2 . This means that each of ^ and T™a N (n = 1, 2, . . . , 2 Ar_2 ) is a positive matrix, and 

hence Eq. (|24| indeed takes the form (fTB")) . 

If there are degeneracies in some of the eigenspaces of the matrices p]^ edia and p^odia' there is a possibility that we 
cannot choose a common state \n) that represents the simultaneous eigenstate \n, Pi,vn) for p, v = 0,x,y,z (Fig. |4j. 
Let us then inspect the degeneracies in more detail. 

Suppose that the matrices pjf edia and p^odia snare the eigenstate \no, pi, vn) with the respective eigenvalues A™ 
and r™°. We can choose the state | no, P\,vn) even when each of the eigenvalues A™ and r™° has degeneracies 

in its own eigenspace. Suppose also that the matrices p^ edia and p^ dia share the eigenstate \n\, p[, vn) with the 
respective eigenvalues A™, 1 and t" 1 . After close inspection, we can state the following: if the eigenvalues r"° and r" 1 
are not degenerate, the states |no, pi,^/v) and \n\, p' 11 i>n) are orthogonal to each other. The only possibility that 
we cannot choose a common state \n) then occurs when all the matrices p^ edia (p = 0, x, y, z) have degeneracies in 
the corresponding eigenspaces and/or all the matrices p^edm ( v = 0,x, y, z) have degeneracies in the corresponding 
eigenspaces. 

We can thereby break down the whole eigenspace into blocks. We form a block of eigenspace in which all the matrices 
/°media = 0, ac, y, z) have degeneracies (Case A) and/or all the matrices p^odia ( p = ^i x iV^ z ) have degeneracies 
(Case B). Let us denote each block as T-Ln with the dimensionality Dn- Let us choose an arbitrary orthonormal set of 
states \n)n (n = 1, 2, . . . , D^) in the block Ha- Then we sum the terms ^2 A"cr^ ® S„ T "°1v inside each block Ha to 
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have 



E E A X® E t >n 

■ \rLj n&'Hn H=0,x,y,z u=0,x,y,z 



tr 1JV (e-^( CT 0|i)-(i| fi ) g,tr 1JV e-^(^)E„=i |n)«(n| fi in Case A, 

( n V ( 2?) 

triiv e-^(^) E^i |»>ft<n| fi ® triw(e-^^)|l) ft <l| ft ) in Case B. 



This shows that Eq. (|24p can be summarized into the form Eq. (|16[) . where the summation in the right-hand side of 
Eq. ([To]) is taken over the blocks n. Thus, Theorem 1 is proved. □ 

The contraposition of Theorem 1 gives us a necessary condition to generate the entanglement between the focused 
spins by the interaction ©. We can see from the proof that we can extend this theorem to be applicable to any 
bipartite systems which indirectly interact with each other at arbitrary temperatures. Let us consider the interaction 
between the bipartite system with C levels; namely 

£ 2 -l 

Hint = E (^1 ® ^Lnodia + ® #W>edia) + ^mcdia, (28) 
i=l 

where {a\, a^jf^ 1 are the bases of the bipartite systems and we define {a®, cr^-} as the identity matrices. In Eq. (|18l) . 
we had //, v = 0, x, y, z for the two s = 1/2 spins, while we have (i, v = 0, 1, 2, . . . ,L — 1 for this bipartite system. 
The extension of the proof is straightforward. 



B. Effective fields 



We next show that the effective fields on the focused spins have serious effects on the generation of the long-distance 
entanglement. In order to make the discussion clear, we first give the formal expression of the effective Hamiltonian 
of the spins 1 and N. 

We investigate the ground state of the total Hamiltonian iJ to t of the system in Fig. [5] (a). For this purpose, we 
separate the Hamiltonian ([]} as follows: 

Htot = Hq + Hi, Hq = -ffmedia, 

Hi = -ff couple + Hi,p. (29) 

We define the ground state of Hq as |?/>" lcdla ) with the eigenvalue £ , ™ edla . Because we assumed that the ground state 
of the media system is not degenerate, we also define the first excitation energy of iT me dia as <5i?f lcdla (> 0). 
We then assume 

\\Hi\\ <<5^5 nedia (30) 
and consider the term Hi as perturbation. The unperturbed ground state is given by 

P = /lJV® W° dia )W Cdia |: (31) 

where Juv is the identity matrix of the spins 1 and N. 

In order to define the effective Hamiltonian, we consider a Green's function 

F p H P = P° F l H Po (32) 

& — -TO-ntot-n) & — H tot 

and define the effective Hamiltonian as 

(33) 



E - P Q H tot P E-H° s ' 
We calculate the effective Hamiltonian by expanding Eq. ([52)) with respect to Hi to have 

P ° p TT~ P ° 

Hj — -Htot 

= Pq E-H -Hi P ° 
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Wc can rewrite Eq. (|34|) as 

P °E - H tot P ° = E - P H P - PoH.Po - PqH^-^QqH^ + °^ H ^ )' (35) 

where Qo = /tot — Po with 7 t ot the identity matrix of the total system. We can indeed confirm Eq. Q35p by expanding 
it. The pole of the Green's function §5§ is given by E = Eff edia + 0(||.Hi||), and hence we substitute E = £™ odia 
into the term P HiQq e-h QoHiPo- We thereby obtain the effective Hamiltonian as 

H cff = P H Q P Q + PoH^o + PoH.Qo — ^ - QoPiPo (36) 

Let us calculate the effective Hamiltonian (|36|) using (|29|) . Because the Hamiltonian Plf and the projection 
operators Po and Qo commute with each other, we have 

P PiP = H LF ® IC° dia ) (C cdia | + PoH coupic P (37) 

and 

1 1 

PqHxQo gmedia. _ Hq QoH X Pq = P P C ouplcQo dia _ — Qo^couplc^O, (38) 

where we utilized the equation PqQo = Po(^tot — Po) = 0. Then, we obtain the effective Hamiltonian for the spins 1 
and N, which we refer to as Hf%, up to the constant P™ edla , in the form 

Ht% = (Vr dia |^ eff |C edia ) = H LF + (^ 10dia |^couple|V'o nCdia > 

+ (V^^'Igcoup^Qo Fmcdi ^ „ Qogcoupk# mCdia )- (39) 

Note that the effective Hamiltonian can be divided into the part which includes only Plf and the part which includes 
only P^ C oupie- The expression (|3^|) is essentially the same as the one derived in Ref. @. 

In order to relate the effective Hamiltonian (|39[) to the argument in Sec. UH we prove the following equation 

tr lw J-e-^ = J- c -^ +0(||P 1 || 2 ) (40) 

in the low-temperature limit j3 — > oo, where tr^ denotes the trace operation on the system except the focused spins 1 
and TV. Equation (|4"0")) means that the effective Hamiltonian defined in Eq. (|3"9"|) gives the density matrix up to the 
second order. The proof is given as follows. First, Eq. (j3"3")l is Laplace-transformed to the following equation: 

Because of the inequality f[5U|) . perturbation theory yields 

-J-e-? 11 ™ = -!-(Po + Qo)e~ m ^(Po + Qo) 

= -!-Poe-' 3H ^Po + O(\\H 1 \\)(P A 1 Q + Q A\P ) + 0(||I?i|| 2 ), (42) 

^tot 

in the low-temperature limit ft — > oo, where we define A\ as an O(l) operator of the total system. We can obtain 
Eq. (gnj from Eqs. (gTJ and |42]). 

Let us characterize each term of the effective Hamiltonian in Eq. (f3"9"]). The first term in Eq. (|3"9")l is the local 
fields on the focused spins, whereas the second term gives the effective fields of the order of ||Pi||; the second term 
does not generate the interaction term because # COU pie given in Eq. flSJ does not include the terms such as a\ £§> a 3 N 
= x,y,z). The third term gives the effective interaction between the focused spins as well as the effective fields 
of order ||Pi|| 2 . The effective Hamiltonian is thereby summerized as follows: 

Hf N = H LF + J2 (^1+^^)+ E JlM®^, (43) 

i=x,y,z i,j—x,y,z 
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where 

E (hT s 4 + /4 eff <4) = (C odia l^cou Pl oic cdia ) + odmf), 

i—x,y,z 

E J>i®4t = 0(\\H 1 \\ 2 ). (44) 

We show that a necessary condition for the media system to generate the long-distance entanglement is that the 
fields of the order of \\Hi\\ vanish, namely, 

(Vo mcdia |ffcoupic|^o n ° dia ) + Hlf = 0. (45) 

Note that the condition (|45|) does not ensure the maximum long-distance entanglement, namely C(pin) = 1; the 
effective fields of the order of ||-ffi|| 2 can still exist even if the condition (j45|) is satisfied, which are of the same order 
as the indirect interaction in the limit of ([7]) . When we explicitly know the form of the effective Hamiltonian, we can 
cancel the effective fields by choosing the local fields on the spins 1 and TV as 

H LF = E (h\ai+h%a j N ) = -W cdia |i?coupie|^o m ° dia >- (46) 

i— x,y ,z 

Let us consider the following special case: 

#couplc= E ( J A a \® a3 2 + J 'B 3cJ N®V j N-l) (47) 

i,j=x,y,z 

with the condition 

(^o modia |^coupic|^ lodia ) = 0. (48) 
In this case, we do not need i?LF to satisfy Eq. (|4"S"|) . namely we can put 7?lf = 0. The condition (|48[) is reduced to 

E W Cdia |^^2l^o ncdia )-0 and E W° dia |^'<-il^o ncdia >=0 (49) 

j=x,y,z j=x,y,z 

for i = x, y, z. This occurs in any systems with the time-reversal symmetry; quantum spin chains without external 
fields and other odd-body interactions typically satisfy Eq. (|49|) , even if the Dzyaloshinskii-Moriya interaction (29l.l30j 
is included in the spin-spin interaction. 

In the derivation (f59"| of the effective Hamiltonian, we assumed that the effective Hamiltonian can be obtained by 
the second-order perturbation. However, there are cases in which we cannot obtain the effective interaction by the 
second-order perturbation. For example, in the cases where the indirect interaction is 'classical,' we cannot obtain the 
effective interaction by the second order; if the effective Hamiltonian for a classical interaction would be given as in 
Eq. (|39p . we could always obtain the maximum entanglement by choosing the local fields {h\, h l N }i =x . y . z properly [25| . 
but this is contradictory to Theorem 1. If the second-order perturbation vanishes, the effective interaction may depend 
on the local fields {h\, h l N }i =x . y . z and takes a more complicated form in higher-order perturbations. 

IV. THE GENERATION OF THE LONG-DISTANCE ENTANGLEMENT BY THE LOCAL FIELDS 

In the present section, we discuss the generation of the long-distance entanglement by the use of only the local 
fields. Throughout the present section, we consider the system ^ in Fig. [21(b). 

A. Effective Hamiltonian 

We can realize the condition equivalent to (|7|) in the system ([§]) by increasing the amplitudes the local fields of 
{h l 2\i =x ,y,z and {h* N _ 1 }i =x ^ y ^ z with the local fields {h\}i =x ^ y ^ z and {h? N }i =x;y ^ z canceling the resulting effective fields. 

In the limits \h2\ — > oo and |/i_/v-i| — > oo, the total Hamiltonian ((SJ) can be transformed into the following effective 
Hamiltonian after tracing out the spins 2 and N — 1 : 

rreff _ \ ' / ji,j,eS i j , ji,j,eS % j \ 
n tot — 2^, V J A a l <7 3 + J B <7 N-2 a N) 

i,j=x,y,z 

+ E + + + ^ (50) 

i— x,y,z 
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equivalent? 





spin 1 spin 3 



FIG. 5: Schematic picture on the equivalence between the decrease of the interaction and the increase of the local fields. We 
can effectively weaken the indirect interaction between the spins 1 and 3 by increasing the local field on the spin 2. 



The coupling parameters {J l A' cS }i,j= x ,y,z and { J l g' cS }i,j= x . y , z in Eq. ([50]) approach zero as | /i2 1 1 and | Hiv— i | — 1 • 
The amplitudes of the effective fields {h\ e }i =x , y ,z and {h l £ }i =x , y ,z are also of the order of 1 7^2 1 1 and \Kn-i\~ 1 , 
which is the same as the coupling parameters {J^' e }i,j=x,y,z and {Jg*' e }i,j=x,y,z- We therefore apply weak local 
fields {h l { eS }i=x,y,z and {h l ^ S }i =x ,y, z of the order of | 1 1 and \Jin-i\~ 1 - In order to achieve the condition (|4"5*j) . we 
apply the local fields hi = —hf and 1%n = —h c ^ on the spins 1 and N to cancel the effective fields, so that 



i—x,y,z 

E (/4 + 4 e£I Kr = o. 



(51) 



In the following, for simplicity, we let the coupling parameters {J]\ ', J % a , Jb" j JB J }i,j=x,y,z ^ e ^ ne XYZ interaction 
and {h\, h? N _- i }i= x ,y, z applied only in the z direction: 

tftot = E ( J A (T l <T 2 + J\ a \ a \ + ^B fT 7V-2 f7 ]v-l + J B IJ N-1< T 

,y,z 

E (hi a i + h l N a l N ) + h\a\ + h z N _ 1 a z N _ 1 + H medi 



N) 



i—x,y,z 



(52) 



-x,y,z 



In Section IIVB1 we show the decrease of the effective coupling interaction mathematically and discuss how precisely 
we can achieve the condition ([7} by increasing [ /i-a | and | /ijv— 1 1 • 



B. The equivalence between the decrease of the interaction and the increase of the local fields 

We here show the equivalence between the decrease of the interaction and the increase of the local fields. In order 
to study the mathematical structure generally, we consider the quantum system shown in Fig. [5] instead of the system 
in Fig. [5] (b); the spins 1 and 3 indirectly interact with each other through the spin 2, while the spin 3 is coupled to 
the environmental system. We define the environment as an arbitrary n en v-dimcnsional quantum system; that is, the 
dimensionality of the total system is 8n en v The total Hamiltonian is given as follows: 

,cnv ? 

H%f* = E (JVK + J'vi4), (53) 

z— x,y,z 

where H^^ cnv is an arbitrary 2n en v-dimensional Hamiltonian of the spin 3 and the environment. The effective Hamil- 
tonian Hf s 3 env after tracing out the spin 2 takes a simple form in the limit |/if | — > oo. We apply the formula (|36|) to 
Hamiltonian (|53[) to obtain the form of Hf 3 onv . 

Theorem 2: The effective Hamiltonian cnv for the Hamiltonian (|53p is given as follows in the limit of — > 00: 

^i cff 3,cnv= E J^^l+H^ + hl^al + h^al (54) 
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where 



jx j x jy j y 

jx,eS _ J J n/.eff _ J J 7Z,cff _ n 



jx jy j x 7 y 

hl' cS = -J z - h z 3 ' cS = -j' z - (55) 



In addition, the Hamiltonian Hf s 3 cnv satisfies 



tr 2 ^ + Od^r 2 ) (56) 

^tot -^tot 

in the limit of f3 — > oo, where tr 2 denotes the trace operation only on the spin 2. 

Comment: The indirect interactions along the x and y-axes decay as |/j.f| in the limit \h 2 \ — > oo, while the indirect 
interaction in the z-axis decays as |ft|| -t with t > 2. In other words, the indirect interaction in the direction of the 
local field on the spin 2 decays more rapidly. 

Proof: In order to prove the present theorem, we follow the same calculations as in Section UlIBI We separate the 
Hamiltonian (|53p as follows: 

-fftot = Ho + H±, H = h 2 a 2 , 

Hi = i^i3. l 2 Pl ° + H 3yCnv . (57) 

We define the ground state of H as ^2) with the eigenvalue ~h 2 . 

First, we prove Eq. (fMf. We can calculate the effective Hamiltonian as in Eq. (|36p : 

H cS = P H Po + PoHiP + PoffiQo , z 1 „ QoHiPo, (58) 

— n 2 - -Ho 

where Pq = |t2)(t2l®^i,3,env with Ii^^nv the identity matrix of the spins 1, 3 and the environment, and Qo = hot — Pa- 
Because H 3 ^ onv commutes with Pq, we have 

H cS =tt 2 )(U® 

-K + H h3 , cnv + (t 2 | (H£f° + H™t c Qo _ h ^_ H Q»H?*t) |ta>] • (59) 

The effective Hamiltonian H° s 3 onv of the spins 1, 3 and the environment is thereby given by H° B 3 cnv = (T2 |-ff cff IT2) ■ 
By calculating Eq. (|59|) up to the constant component — h%, we arrive at Eq. ([54)) with Eq. (|55|) . 
Next, we prove Eq. ([55)1 . Assuming h 2 ^> \\Hi\\, we obtain 

p— P — /3-fftot 

-s = ^0-^ Po + 0(|^r 1 )(P AiQ + Qo4^o) + 0(|^r 2 ) (60) 

in the low-temperature limit as in Eq. (112"1) . where we define as an O(l) operator of the total system. We then 
obtain Eq. (|56j) from Eqs. (|4Tj) and (|60l) . This completes the proof of Theorem 2. □ 

Applying Theorem 2 to the Hamiltonian ((9]), we obtain the effective Hamiltonian H^ t of the form ([50jl . For the 
Hamiltonian (|52p . more specifically, the effective Hamiltonian H°^ t is given by 

Htot =HhF + -^LF + ^couple + Pmcdia + ^nicdiai (61) 



where 



eff _ I jz ,J A' J A \z I ~rz I 



jx j'x TV jV jx j'x jV jV 

rreff J A J A „x x J A A „y Jy J B J B x x J B B „y „y 

-"couple — ZI~ (T 1 (7 3 J- z a l a 3 Ti <J N-2 a N 7~z a N-2 a N' 

n 2 <i 2 ,l N—l n N—l 

H e jL dia ^-(ji+^yi- (j B +^y N - 2 . (62) 

We thus achieve ||-ff C oupic| in the limits \h 2 \ — > oo and | /i^r 1 1 —> oo. The local fields h\ and h z N which cancel the 

effective fields arc given by 

jx jy jx jy 

hl = Jl + -4^, h z N = r B+ J -pl. (63) 
n 2 n N-i 
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(a) (b) 

FIG. 6: We plot (a) the long-distance entanglement C(puv) and (b) the first excitation energy Sei against the field ho in the 
three cases 7 = (solid line), 7 = 0.03 (broken line) and 7 = 0.05 (chained line). 
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0.03 


0.05 


Hamiltonian 


J64j) 


0.9832 


0.6743 


0.1631 


Hamiltonian 


§B 


0.9834 


0.6745 


0.1632 



TABLE I: Comparison of the entanglement for ho = 100. The effective Hamiltonian of the total Hamiltonian (|64p is equal to 
the Hamltonian ((65} by the second order perturbation. The difference between them is of the order of h^ 2 = 10~ 4 . 



C. Numerical demonstration 



In the present subsection, we numerically demonstrate the generation of the long-distance entanglement by the use 
of the local fields for an XY spin chain. We consider the total Hamiltonian 

Htot =(1 + 7)(< J l '2 + °2°3 + °N-2 a N-l + VN-iVn) 

+ (1 - 7)Ol02 + °t°% + ^N-2 a N-l + a N-l aV N) 
(1 -7 2 ) 

+ j K + <J Z N ) + h {<T Z 2 + 0-^_ x ) + ff m edia, 



ho 



N-3 



H media = [(1 + 7)of of +1 + (1 - 7)^Vf +1 ] + ^-fr^-(^ + ^-2), (64) 

ho 



1=3 



where we choose the local fields on the spins 2 and N — 1 as h\ = h z N _ 1 = ho- We introduced the third term in H to t 
accordingly to Eq. (|55|) in order to cancel the effective fields on the spins 1 and N which are generated after tracing 
out the spins 2 and N — 1. We also introduced the second term in 7? me dia to cancel the effective fields on the media 
spins 3 and N — 2 for simplicity. We obtain the effective Hamiltonian H^ t of the form (|5T|) as 

JT off — (I+7) . . (1 — 7) f y y y y s 

iV-3 

+ E [(! + °f+i + (! - 7)^<i] (65) 

(=3 

in the limit of /io — > 00. The ground state of the XY spin chain without fields satisfies the condition (|48[) because 
the Hamiltonian has the time-reversal symmetry. Note that the indirect interaction through the XY spin chain is 
not a 'classical' one. We can therefore expect that the long-distance entanglement exists in this system in the limit 
of ho — > 00. 

Let us show the numerical calculation of the long-distance entanglement for the system (|64|) with N — 100. In 
Fig. [6] (a), we show the entanglement between the spins 1 and N in the three cases, 7 = 0, 0.03 and 0.05. We can 
see that the entanglement C(p±n) monotonically increases with ho- As the parameter 7 increases, however, we need 
a greater value of ho to generate a large long-distance entanglement. This is because the first excitation energy of 
the media system decreases rapidly as the parameter 7 increases; as is shown in Eq. (|30j). we need to attenuate the 
coupling strength adequately so that it may be much less than the first energy gap of the media system. If the 
condition (|30[) docs not hold, the second term on the right-hand side of Eq. (|40[) may not be ignored and may decrease 
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the purity of the total system, which causes the destruction of the long-distance entanglement. We show in Fig. [6] (b) 
the first excitation energy Sei of the total system, which decreases as h^ 2 in the limit of h a — > oo. This would make 
the long-distance entanglement fragile to thermal fluctuation. 

The value of the long-distance entanglement is almost the same if we consider the total Hamiltonian (|65p with 
N = 98 from the beginning. We compare the values in Table. U for ho = 100. The difference of C(pin) is of the order 
of hg 2 as is expected from Eq. ([56]) . 

V. FINITE COUPLING 

In the previous sections, we consider the weak coupling limit (J7J in order to generate the long-distance entanglement. 
By making the coupling strength to zero, however, the first excitation energy of the total system vanishes and the 
long-distance entanglement becomes extremely fragile against thermal fluctuation. In order to avoid the situation, 
we have try to achieve the long-distance entanglement by as strong coupling as possible. In the present section, we 
consider the Hamiltonian ([T]) in Fig. [2] (a) with a finite coupling Hamiltonian: 

Uncouple 1 1 = const ^ 0. (66) 

We focus on the following two points: 

1. We look for the media system suitable for the generation of the long-distance entanglement. 

2. We do not necessarily need the condition (|45p for the generation of the long-distance entanglement when the 
coupling Hamiltonian is finite. 

First, in Section [V Al we introduce two examples of the XY spin chains with the external fields in the z direction in 
order to demonstrate that the external fields on the media system can enhance the long-distance entanglement with 
the coupling strength fixed. With a finite coupling Hamiltonian, the preparation of a suitable media system is crucial 
for the generation of the long-distance entanglement. For example, in Section HV C[ we showed that the long-distance 
entanglement strongly depends on the parameter 7 if we keep the finite coupling strength. We therefore consider the 
possibilities that we can enhance the capability of the media system with external fields on the media system. 

Second, in Section IV Bl we discuss the case in which a random Hamiltonian, which may generate the effective 
fields (|44[) . is added to the system. There, we tune the coupling strength and see that the coupling strength should be 
neither too strong nor too weak. In the limit of the weak coupling ([7]) of this system, the long-distance entanglement 
always vanishes in any quantum systems which do not satisfy the condition (|45|) . If the coupling has a certain value, 
however, the entanglement may exist in the system without the condition (|45|) . As has been analyzed in Section IlIIBl 
the effective fields on the focused spins are of the first order of the coupling amplitude O(||i? cou pie||)- On the other 
hand, the effective interaction is of the order of ||-ff C oupie|| 2 , and hence the ratio of the effective fields to the effective 
interaction is of the order of ||J? C oupie|| _1 - We therefore expect that by increasing the coupling strength we can 
relatively reduce the effective fields and enhance the long-distance entanglement. 

A. Enhancement of the long-distance entanglement with the external fields 

In the present section, we consider the possibilities that we can enhance the capability of the media system to 
generate the long-distance entanglement by the use of the external fields on the media system. We consider the XY 
spin chains with the external fields in the z direction in the following two cases: random fields and uniform fields. First, 
it has been shown that in the XX spin chain the random fields can enhance the entanglement between short-range 
spin pairs [24[. We thereby expect that a randomness may also enhance the long-distance entanglement. Second, it 
is known that the quantum phase transition occurs in the XY spin chains with the external fields. Then, we expect 
that the long-distance entanglement is enhanced around the critical point because of strong quantum fluctuation. We 
indeed show that the long-distance entanglement is highly enhanced in these two cases. 

We first give the Hamiltonian 

fftot =#media + 0.02(afaf + afa y 2 + a x N _ Y a x N + (7 v N _ l( j v N ), (67) 

where 

N-2 N-l 

ffmcdia - £ [a(l + 7 K*f+i + (1 - 7)*K + i] + E h '°t (68) 

1=2 1=2 
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FIG. 7: The distribution of (a) the long-distance entanglement and (b) the first excitation energy. We determined each of the 
external fields {/i;}fl 2 stochastically out of the uniform distribution [—1.5, 1.5], to obtain 10 7 samples. In (a), the value for the 
bin from to 0.002 is 0.134, which is out of the range of the plot. The average of the entanglement is 0.162. Without random 
fields, the value of the long-distance entanglement would be 1.69 x 10 -4 and the first excitation energy would be 1.19 x 10~ 8 . 
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FIG. 8: The entanglement against ho, (a) for 7 = 0.3 (solid line), 7 = 0.5 (broken line) and 7 = 0.7 (chained line) with N — 100 
and (b) for N = 50 (solid line), N = 100 (broken line) and N = 200 (chained line) with 7 = 0.5. We apply the uniform fields ho 
on the XY spin chain (|67|) . The long-distance entanglement can be enhanced around the point ho = 27. As the spin number 
TV increases, the peak approaches the critical point ho = 27 and becomes sharp. 



and a is an integer which has the values ±1. The property of the media system qualitatively changes depending on 
the sign a as is shown in the following. In the present system, the indirect interaction between the spins 1 and N 
depends on the external fields {hf}]^ 1 . We then expect the possibilities that {hf }j^J 2 1 enhance the generation of the 
long-distance entanglement. 

First, we apply the random fields {hf distributed uniformly in the range —1.5 < hf < 1.5 (I = 1, 2, . . . , N— 1). 
We let a = +1, 7 = 0.1 and TV = 100. For hf = (I = 2, . . . , 99), the long-distance entanglement C(pin) and the 
first excitation energy 8e\ are equal to 1.69 x 10 -4 and 1.19 x 10~ 8 , respectively. In Fig. [TJ we show the distributions 
of C(pin) and log 10 Sei in the case of the random fields. The average of C(puv) is 0.162; we can see that the 
random fields highly improve the entanglement and the first excitation energy. There are even the cases in which the 
entanglement C(puv) is more than 0.95; the largest value of the entanglement is 0.982. This suggests that we can 
significantly enhance the long-distance entanglement by choosing the external fields optimally. 

Second, we apply the uniform fields in the case a = — 1: 

hf = h (69) 

for I = 2, . . . , N — 1. In this case, a quantum phase transition occurs at the point ho = 2j in the limit of N — > 00 [3l"l ]. 
In Fig. [8] (a), we show the plots of the long-distance entanglement against the parameter ho for 7 = 0.3, 0.5, 0.7 and 
TV = 100. We can see that the entanglement is highly enhanced near the critical point ho = 2-f. In Fig. [5] (b), we 
show for 7 = 0.5 that the peak becomes sharp and close to the point ho = 27 as the spin number N increases from 
50 to 200. For a = 1, we cannot achieve such enhancement by the external fields even at the critical point. 
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FIG. 9: The distribution of the entanglement C(pijv) (a) for J = 0.05, (b) for J = 0.25 and (c) for J = 0.5 for the Hamilto- 
nian (|70l) . For the panels (a)— (c), we used 250 000 samples for each case. In the panel (d), we plot the average of C(puv) for 
each value of J with 2500 samples for each point to calculate the average of the entanglement. 



J 


0.05 


0.25 


0.5 


Entanglement without random Hamiltonian 


0.973 


0.551 


0.102 


Entanglement with random Hamiltonian 


0.03 


0.27 


0.07 



TABLE II: The average values of the entanglement without and with the random Hamiltonian in Eq. (|70[) 



B. Small spin chains with the random Hamiltonian 

Next, we consider a quantum system which docs not satisfy the condition (|45|) . In order to discuss the effect of the 
breaking of the condition (|45[) . wc consider the following spin chain with N = 8, 

H = £>f af +1 + of af +1 ) + J(aM + °X°l + « + «) 

1=2 

+E E <"T^H-i. ( 70 ) 

1=1 0,x,y,js 

where er° denotes the identity operator. We add the random coupling {<5J; IJ }i J= o,a;,y,z (I = 1,2,..., 7) to the XX 
spin chains; we assume that the randomness is uniform in the range —0.05 < <5J/ J < 0.05 for i,j = 0,x,y,z and 
/ = 1,2, ... ,7. In this system, the condition (|4"5")) is not satisfied, and hence the entanglement C(pxjv) is equal to zero 
in the limit J 0. 

In Fig. [S] (a)-(c), we show the distribution of the entanglement C(pijv) for J = 0.05, 0.25 and 0.5. In Table HT1 
we show the averages of the entanglement in these three cases. We also show the values of the entanglement without 
the random Hamiltonian, namely J\' 3 = for i,j = 0,x,y,z and I = 1,2,..., 7. From Table ITT1 we can see that 
the destruction of the entanglement by the random Hamiltonian is smaller as the coupling parameter J increases. 
However, a too strong coupling parameter does not give the optimum value of the entanglement nor does a too weak 
coupling parameter. If the coupling parameter J is small, the destruction of the entanglement by the effective fields 
is serious. On the other hand, if the coupling parameter J is large, the entanglement between the focused spins and 
the media system becomes non-negligible, the purity of the focused spins decreases, and hence the entanglement is 
destroyed. The entanglement destruction is thus different physically in these two cases. These two effects compete in 
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giving the optimum value of the coupling parameter. As is shown in Fig. [9] (d), we can achieve the maximum average 
of the entanglement with the coupling parameter J ~ 0.25. 

VI. SUMMARY AND CONCLUSION 

We have analytically and numerically studied the generation of the long-distance entanglement by the use of the weak 
coupling ([7]). We gave the two necessary conditions for the media system to generate the long-distance entanglement. 
The first one is that the indirect interaction must not be 'classical' as has been defined in Eq. (fTDj) . The second one is 
that the effective fields of the first order have to vanish as in Eq. (|45|) . The first condition, in particular, is applicable 
to any bipartite systems indirectly interacting with each other at any temperatures. The second condition can be 
satisfied artificially by applying the local fields so as to cancel the effective fields. In this sense, we may overcome 
the constraint of the second condition if we know the effective fields on the focused spins. We have shown that these 
two conditions can be satisfied in various systems; this means that many quantum systems have potentials of the 
generation of the long-distance entanglement. 

Next, we discussed the generation of the long-distance entanglement by the use of only the external fields. As 
shown in Fig. [2] (b), we have to control the local fields on the focused spins and the part of the system adjacent to 
the focused spins. We have to apply strong fields on the adjacent spins and the weak fields on the focused spins. 
Then, we achieved the condition mathematically equivalent to Eq. (|7|). The strong fields contribute to attenuate the 
coupling between the focused spins and the media system and weak fields is necessary to cancel the effective fields on 
the focused spins. Because we utilize the fields with finite amplitudes, we can achieve the condition ([7]) by the second 
order approximation as shown in Eq. (|56|). but the degree of accuracy rapidly improves as the amplitude of the strong 
fields increases. Our result also makes it possible to control the interaction parameter by the local fields, and hence 
it may be also applicable to accurate control of the interaction. 

Finally, we showed the cases where the coupling strength is non-zero. We first discussed the two cases in which the 
external fields on the media system can enhance the generation of the long-distance entanglement; we proposed the XY 
spin chain with the uniform fields and the random fields with the Hamiltonian (|67[) . The long-distance entanglement 
is enhanced by the random fields on average; there are a few cases in which the entanglement is enhanced to be nearly 
equal to 1. This means that we can improve the capability of the media system to a great extent. The long-distance 
entanglement is also highly enhanced near the point ho = 27, which corresponds to the point of the quantum phase 
transition in the limit N — > 00. Second, we discuss the cases in which the effective fields of the first order remain to be 
left; we consider the XX spin chain with random interactions. In this system, the entanglement always vanishes in the 
limit of the weak coupling (|45[) . A too strong coupling causes the decrease of the purity because of the entanglement 
with the media system, while a too weak coupling increases the relative amplitude of the effective fields to the effective 
interaction. The optimal coupling strength is determined so as to make the both effects minimum. If we choose the 
coupling strength properly, the long-distance entanglement still remains to some extent. 

In conclusion, we have generally researched necessary conditions for the generation of the entanglement. Our results 
show that the conditions for the generation are not so strict in the ground states. For the practical application, however, 
we cannot consider the low-temperature limit because it is not realized experimentally, and hence the amplitude of 
the first excitation energy is important. We should keep a certain coupling strength in order to increase the first 
excitation energy as much as possible; in the limit of weak coupling, the first excitation energy always vanishes. For 
this purpose, it is essential to find or construct a media system best suited to the generation of the long-distance 
entanglement. We present the possibilities that the capability of the media system can be improved by the external 
manipulation. However, we do not understand the principle on how to improve the capability of the media system. 
Then, the analysis of the capability of the media system will be the next problem for the practical application of 
the long-distance entanglement. There are several trials [l(| EH aimed at finding efficient media systems for the 
entanglement generation. 
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Appendix A: Relation between Theorem 1 and the quantum discord 

In Theorem 1, we show a necessary condition for the entanglement to be generated via a media system. In the 
present section, we answer the following question: if the condition (|10[) is satisfied, can the quantum discord still 
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exist? The answer to this question is yes and we show an example in the following. 

First, we review the definition of the quantum discord (32[. Let us consider a two spin system with the density 
matrix puv- The quantum discord Q(/9ijv) between the spins 1 and N is defined as follows: 

Q(pm)=l(pi N )-J(p 1N ), (Al) 

where I{p) is the quantum mutual information defined by 

l(p) = S( Pl ) + S(p N ) - S(p 1N ) (A2) 

with S(p) the von Neumann entropy S(p) = tr(plnp). On the other hand, J(p) is the optimized classical mutual 
information, which is the maximum information obtained from the measurement of the spins 1 or N, and is defined 
by 

J(p) = S(p N ) - mm^pjSipN^), (A3) 

3 

where S(pn) is the initial von Neumann entropy of the spin N and y"V PjS{pM\Yij ) is the average of the von Neumann 
entropy after the measurement of the spin 1 in the basis of ILy. If the quantum discord (|A1|) has a non-zero value, 
the correlation between these two spins may not be explained by classical theory. 
Let us consider the Hamiltonian 

H tot = H int + H hF , (A4) 

where 

-Hint = + cr 2 cr 3> 

H LF = J2 (h\4 + hi4)- (A5) 

i— x,y,z 

The Hamiltonian (IA5I) satisfies the condition ffTOl) as 



[H A (a 1 ),H B (a 3 )}=0, 

H A = ofo%, H B = a%al (A6) 

Therefore, the entanglement can never exist between the spins 1 and 3 in the ground state of -f/tot for any values of 
the local fields {h\, h l 3 }i =x ^ y ^ z . Indeed, the density matrix P13 is given by 



pi 3 = lim tr 



/3-T^ 13 VZ tot (/?) 



e -/3ff,ot 



/0.7286 0.1250' 

0.1250 0.1250 

0.1250 0.1250 

V0.1250 0.02145y 



(A7) 



for h\ =h\ = \ and hf = = h\ = h\ = 0. This system has no entanglement. 

However, it has a non-zero quantum discord. We utilize the criterion in Ref. (33j to prove this. First, we separate 
the density matrix into the following four blocks: 

n _ A).7286 \ 12 _ / 0.1250\ 
9 ~ \ 0.1250j ,P ~ ^0.1250 )' 

21= f 0.1250\ 22 _ /0.1250 \ 
p ^0.1250 J ,p \ 0.02145/ ' y ' 

A necessary and sufficient condition for the zero discord is given by the following two statements: 

[p« (p«)t]=0fori,j = l,2 (A9) 

and 

[p ij , /^']=0 fori, j,i',j' = 1,2. (AlO) 



19 



The density matrix (|A7|) satisfies the first condition (| A9|) because it is a real matrix. However, the second (|A10[) 
condition is not satisfied. Indeed, 

11 12 _ ( 0.09107\ i2i\_( 0.01562\ , 
P P - ^0.01562 ) ' p p ~ ^0.09107 ) ' [AU) 

and we have p 11 p 12 ^ p 12 p u . Therefore, there exists a quantum discord between the spins 1 and 3. This shows that 
the condition in Theorem 1 is applicable only to the existence of the entanglement, not the quantum discord. So far, 
we are not sure whether there exists a condition for the indirect interaction to generate a quantum discord. 
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